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In preparation for future, large-scale, multi-object, high-resolution spectroscopic surveys of the Galaxy, we present a 
series of tests of the precision in radial velocity and chemical abundances that any such project can achieve at a 4m 
class telescope. We briefly discuss a number of science cases that aim at studying the chemo-dynamical history of the 
major Galactic components (bulge, thin and thick disks, and halo) - either as a follow-up to the Gaia mission or on their 
own merits. Based on a large grid of synthetic spectra that cover the full range in stellar parameters of typical survey 
targets, we devise an optimal wavelength range and argue for a moderately high-resolution spectrograph. As a result, the 
kinematic precision is not limited by any of these factors, but will practically only suffer from systematic effects, easily 
reaching uncertainties <lkms . Under realistic survey conditions (namely, considering stars brighter than r = 16mag 
with reasonable exposure times) we prefer an ideal resolving power of R ~20000 on average, for an overall wavelength 
range (with a common two-arm spectrograph design) of [395;456.5] nm and [587;673] nm. We show for the first time on 
a general basis that it is possible to measure chemical abundance ratios to better than 0. 1 dex for many species (Fe, Mg, 
Si, Ca, Ti, Na, Al, V, Cr, Mn, Co, Ni, Y, Ba, Nd, Eu) and to an accuracy of about 0.2 dex for other species such as Zr, La, 
and Sr. While our feasibility study was explicitly carried out for the 4MOST facility, the results can be readily applied to 
CNj ■ and used for any other conceptual design study for high-resolution spectrographs. 
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The formation and evolution of a galaxy like the Milky Way ers the Ca 11 tv[ plet region, 847-874 nm (|Katz et al.||200i 



(MW) still poses many puzzles to modern astrophysics. The so tha t only for a fraction of the brighter stars, limited in- 

Gaia satellite, expected to be launched in 201 3, will make an formation about the chemical composition will be obtained, 

important step towards addressing many of those questions Moreover, it will not be possible to measure detailed ele- 

pertaining to Galactic evolution (Munari & Castelli 2000; ment abundance ratios for a wide range of elements with 

Munari 2001 ; Perryman et al. 2001 ; Katz et al. 2004; Bailer- Gaia - 

Jones 2002; Lindegren 2010). In the first place, Gaia will while Gaia will provide a three-dimensional map of up 

deliver accurate stellar distances from parallaxes and, cou- t o a billion stars in the MW (about 1% of the total num- 

pled with proper motion measurements, transverse veloci- ber of stars) and thus revolutionise the understanding of our 

ties for more than 10 9 stars as faint as V~20mag. Radial Galaxy, it is not sufficient in itself to address a number of 

velocities can also be obtained from Gaia's spectrograph, important questions in Galactic evolution. 

although this will be limited to approximately 1.5 x 10 8 _ . .,, , _ , , . , ,. . 

Gaia will benefit hugely from additional, ground-based 

. " \ ., _ , spectroscopic information (detailed chemical abundances 
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p e l low " ' and precise hne-of-sight velocities) for a broad range of 

stellar types, thus probing the various populations of the 
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MW. Large survey missions that are currently designed as 
follow-ups to Gaia, where it lacks spectroscopic sensitivity, 
come with their own merits and allow us to gain ever deeper 
knowledge about the formation and chemo-dynamical evo- 
lution of the MW. 

Present day multi-fibre facilities^ such as the European 
Southern Observatory/Very Large Telescope (ESO/VLT) 
multi-object facility FLAMES (140 fibres; Pasquini et 
al. 2002), Australian Astronomical Observatory's (AAO) 
AAOmega+2dF instrument (—400 fibres; Sharp et al. 2006), 
Hectoechelle (300 fibres; Szentgyorgyi et al. 2011) and 
Hectospec (300 fibres; Fabricant et al. 2005) on the MMT, 
MIKE-fibers (128 "red" and 128 "blue" fibres; Bernstein et 
al. 2003) on the Magellan telescope, or the 6dF instrument 
(150 fibres, Watson et al. 1998) on the Schmidt telescope of 
the AAO, currently engaged in the RAdial Velocity Experi- 
ment (see Steinmetz et al. 2006), can provide large numbers 
of spectra and reach intermediate resolutions (R ^20 000 
for FLAMES and R < 8 000 for AAOmega and RAVE), 
but, as of yet, there is no fully operational, dedicated, full- 
time survey of all Galactic components in high-resolution 
mode, in the optical. A first step in this direction is the pub- 
lic Gaia-ESO Survey (Gilmore et al. 2012), which started 
on Dec. 31, 201 1 using FLAMES in high-resolution modes 
to assess the chemo-dynamics of some 10 5 Galactic field 
and cluster stars. The Apache Point Observatory Galactic 
Evolution Experiment, APOGEE at 2.5 m APO, Eisenstein 
et al. (201 1), observes in the near infra-red, in the H band at 
a resolving power of about 30 000. 

In the era of Gaia, a considerable fraction of the as- 
tronomical community is giving thought on how to best 
take advantage of the information flow that Gaia will bring 
along. Much effort is thus devoted to conceive large sur- 
veys to complement and complete the radial velocity and 
chemical pattern information in the MW. For details on 
ongoing surveys and instruments the interested reader is 
referred to RAVE, Steinmetz et al. (2006); Sloan Exten- 
sion for Galactic Understanding and Exploration, SEGUE, 
Yanny et al. (2009); and APOGEE. Details on some stud- 
ies on instrumentation are documented for 4MOST (at the 
4.1m Visible and Infrared Survey Telescope for Astronomy, 
VISTA) in de Jong (2011), de Jong et al. (2012a,b); GYES 
at the 3.6m CFHT, Bonifacio et al. (201 1); High Efficiency 
and Resolution Multi-Element Spectrograph, HERMES, at 
the 4 m Anglo- Australian Telescope, AAT, Barden et al. 
(2010); MOONS, Multi-Object Optical and Near-infrared 
Spectrograph, at the 8.2 m VLT, Cirasuolo et al. (201 1); and 
WEAVE at the 4.2 m William Herschel Telescope, WHT, 
Balcells et al. (2010). In this context, a detailed study of 
the required performance and characteristics of such spec- 
trographs is of prime importance. 

In general, considerations for these fibre-fed spectro- 
scopic survey facilities are driven by a large field of view 



For a comprehensive overview of cur- 

rently commissioned and planned projects see 

|http : / /www . ing . iac . es/weave/moslinks . html| 



to cover large fractions of the sky, by a generous multiplex- 
ing capability and high enough spectral resolving power to 
detect kinematic and chemical substructures throughout all 
major components of the MW (thin and thick disks, bulge, 
and halo), and by a broad wavelength coverage to ensure a 
high measurement precision and the detectability of as large 
a number of elements as possible and to be able to derive the 
stellar parameters and the chemical composition of the stars 
at best. 

What follows refers especially to our effort in un- 
derstanding and defining the characteristics that a typical 
survey-spectrograph should have in order to achieve the pre- 
cision necessary to fulfil typical science cases (see Sect. 2) 
that ground-based Gaia follow-up surveys endeavour to pur- 
sue. While the present tests have been explicitly run as a 
conceptual design study for the 4MOS10 facility, our meth- 
ods and results are generalised and can be readily applied 
for phase A studies of any other, potential spectrograph 
projects in the future. This paper is dedicated to the high- 
resolution mode of the 4MOST survey (and others). An 
analogous study of the low-resolution capabilities of such 
instruments will be presented in a future paper. In particu- 
lar, here we elaborate on our efforts to determine the optimal 
wavelength range and the resolution of any such spectro- 
graph, given the signal-to-noise (S/N) ratios necessary to 
achieve the accuracy requested for typical Galactic science 
cases. 

This paper is organised as follows: in §2 we give an 
overview of typical science cases that drive the design of 
ground-based survey strategies and the requirements these 
pose to a survey instrument. In §3 we discuss the target 
range and spectral types that need to be considered in the 
design studies, while §4 describes the grid of synthetic spec- 
tra we employ for our study. §5 deals with our tests on the 
technical requirements to the spectrograph based on the sci- 
ence questions. We close with the discussion of our design 
study in §6. 

2 General science questions 

The mechanisms of the formation and evolution of the ma- 
jor MW components are encoded in the location, kinemat- 

2 4MOST, the 4 metre Multi-Object Spectroscopic Telescope, is a high- 
multiplex, wide-field (3 degree 2 field-of-view) fibre-fed spectrograph sys- 
tem in its phase A study. It is intended for the ESO/VISTA 4. 1-m telescope 
to work on a dedicated 5-years survey base, stalling in 2019. The spectro- 
graph is designed to provide simultaneously at least 1500 (with a goal of 

3 000) low- (R>5000) and high-resolution (R~20000) spectra. We note 
that the definition of "high" and "low" resolution is often somewhat arbi- 
trary in the literature. To distinguish the two operation modes of 4MOST 
we simply call them high- and low-resolution in the following. In the end, 
it is envisioned to obtain tens of millions of spectra at low resolution and 
more than 1 million spectra at high resolution. The bulk of the high reso- 
lution observations of 4MOST will extend down to an r-band magnitude 
of 16 mag at a signal-to-noise ratio of > 150 A -1 (which is translated in a 
S/N ratio of > 44 pixel -1 at the centre of the blue range with a resolving 
power of 20 000 and a sampling of 2.5) in less than 2-hour exposures. This 
will ultimately allow one to achieve the science goals outlined in Section 2. 
For details on 4MOST see de Jong (201 1); de Jong et al. (2012a, 2012b). 



© 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Wcinhcim 



www.an-journal.org 



Astron. Nachr. / AN (2012) 



3 



ics, and chemistry of their stars: physical conditions of star 
formation suggest that most of the stellar material is pro- 
duced in a finite number of large clumps, characterised by 
some amount of chemical homogeneity. Various dynami- 
cal processes that act upon these stellar aggregates (spiral 
arms, bar, interactions or mergers with satellites) can erase 
the kinematic imprint of objects born together. Also, lighter 
clusters will simply disperse due to the original internal ve- 
locity dispersion and general interaction with the gravita- 
tional field of the MW. 

Detailed stellar chemical compositions can then be used 
as genetic markers to unravel the lineage of common field 
stars. In this context, prime science drivers are inevitably 
the chemical tagging of any chemo-dynamical substruc- 
tures (Freeman & Bland-Hawthorn 2002; Bland-Hawthorn, 
Krumholz & Freeman 2010), including abundance gradi- 
ents and their temporal variations (e.g., Chiappini et al. 
2001), and searches for (extremely) metal-poor stars in the 
halo and bulge such as to characterise the first stellar genera- 
tions in the Universe, thereby constraining galaxy formation 
at the earliest times (e.g., Beers & Christlieb 2005, Karlsson 
etal. 2011). 

For turn-off stars and sub-giants it will be possible to 
derive their ages directly based on the parallaxes measured 
by Gaia, however, for stars on the main-sequence and giant 
stars no age information is available. The age-metallicity 
and age-velocity relations for a stellar population represent 
strong constraints on any chemo-dynamical model. As ages 
are only available for turn-off stars this information will be 
limited to the more nearby stars (about 1 kpc). In practice it 
will still be feasible to determine rough ages using "chem- 
ical clocks" (such as a and neutron capture elements; e.g., 
Tinsley 1976, Matteucci 2003). This can be done with ease 
also at such large distances as the outer disk and the Galactic 
bulge. 

2.1 Thin and thick disks 

Despite recent, major progress in understanding disk for- 
mation in cosmological simulations, many discrepancies re- 
main between the predictions of simulations and the ob- 
served properties of disk galaxies (e.g., Scannapieco et al. 
2009, Piontek & Steinmetz 2011). In fact, mergers may 
not necessarily be the dominant process in the formation of 
disks, but internal evolution processes such as radial migra- 
tion may also play an important role in shaping disk galaxies 
(e.g., Schonrich & Binney 2009, Minchev et al. 2012). 

Recent advances in high resolution spectroscopy have 
shown that a kinematic separation between a thick and a 
thin disk is, in principle, statistically, feasible. The two pop- 
ulations show distinct trends in chemical abundances, but 
when the statistical decomposition is used there is naturally 
an overlap (see, e.g., Fuhrmann 2008; Bensby et al. 2005; 
Reddy et al. 2006; Ruchti et al 2011). It is potentially pos- 
sible to identify the different abundance patterns as the re- 
sult of different star formation histories (Chiappini 2011) 



but they can also be explained as the result of internal redis- 
tribution of the stars in the stellar disk (e.g., Schonrich and 
Binney 2009, cf. Minchev et al. 2012). 

However, an important criticism of current observa- 
tional chemo-dynamical studies of the stellar disk is that the 
samples are confined to the immediate solar neighbourhood 
(e.g., Steinmetz et al. 2006), and that the separation into thin 
and thick disk stars is solely based on the stellar kinematics, 
while radial migration of stars (Minchev & Famaey 2010) or 
accretion (e.g., Williams et al. 201 1 and references therein) 
could bring to the solar neighbourhood stars that were born 
elsewhere contaminating the "local" samples (e.g., Grenon 
1999; Haywood 2008). Stars of a given age would then 
show a spread in metallicity equal to the radial spread in the 
metallicities of the interstellar medium at their time of birth. 
Radial migration effects are thus manifested in the magni- 
tude and evolution of the radial abundance gradients in the 
disk. 

Chemo-dynamical models provide very different pre- 
dictions for several key observables (e.g., age-metallicity 
and age-velocity relations, abundance patterns such as 
[X/Fe] vs. [Fe/H] and the scatter in these relations) at dif- 
ferent galactocentric distances and heights from the plane, 
depending on the details of radial migration and chemical 
evolution model. Moreover, stars from accreted satellites 
could dominate the outer disk, while any population formed 
in-situ will be more important closer to the Galactic plane 
(Villalobos & Helmi 2009). In particular, farther out in the 
disks the predictions are very model-dependent as in those 
regions additional mechanisms are at play, such as star for- 
mation thresholds, accretion, heating and stripping by dark 
matter halos. Therefore, full constraints on the formation 
of this component can only be obtained through joint kine- 
matic and abundance surveys covering the entire Galactic 
disk to a high accuracy. In this context, it has been recently 
demonstrated that the eccentricity distribution of orbits in 
the thick disk can be linked to the accretion vs. heating vs. 
migration scenarios (Sales et al. 2009; Wilson et al. 2010; 
Dierickx et al. 2010). To make significant progress in the 
entire field, it is necessary to map the density of stars in a 
multi-dimensional space, namely comprising age, metallic- 
ity, position and the 3-dimensional velocities. 

2.2 The bulge 

Two main formation scenarios have been proposed to ex- 
plain bulges of spiral galaxies (Kormendy & Kennicutt 
2004). In classical bulges, most stars are formed during 
an early phase of intensive star formation following ei- 
ther the collapse of a proto-galactic cloud or a phase of 
intense merger/accretion events as predicted by cold dark 
matter simulations. On the other hand, boxy or peanut- 
shaped, pseudo-bulges are a consequence of secular evolu- 
tion driven by dynamical instabilities in the inner disk, such 
as the buckling of a bar (Combes et al. 1990). The differ- 
ent scenarios predict different timescales of formation, and 
hence different chemo-dynamical histories for bulges. 
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The formation of the Galactic bulge appears to have 
happened through a combination of processes: high- 
resolution spectroscopic studies of giants al ong low extinc 



tion windows (Zoccali et al. 2006; 2008; IGonzalez et al 



201 lb iHill et al. I l201lb have revealed mostly old 



enhanced stars, compatible with a classical bulge (cf, 
Bensby et al. 201 1). On th e other hand, its boxy shape, re- 



vealed in the near-infrared (Binnev et al. 1997; Dwek et al 



19951) . and its kinematics (e.g., Howard et al. 2008, 2009; 
Shen et al. 2010) suggest that our Galaxy harbours a 
pseudo-bulge. Likewise, the co-existence of two red clump 
populations in particular fields toward the Galactic bulge 
(Mc William & Zoccali 2010) argue for an X-shaped, 
double-funnel shape of the Galactic bulge/bar. 

The power of combining chemistry with kinematic stud- 
ies in the bulge is illustrated by the recent work of Babu- 
siaux et al. (2010), who found two distinct populations in 
the bulge, namely a metal-rich, a-poor population with bar- 
like kinematics, and a metal-poor, a-enhanced population 
that shows rather spheroidal or thick-disk like properties. 
The only way to clarify the critical questions on the current 
structure of the Galactic bulge, its formation and the tran- 
sition to other Galactic components is to map the distribu- 
tion of elemental abundances to high precision and to cou- 
ple these to accurate kinematics over a wide range of bulge 
locations, thereby unveiling the mix of the different popula- 
tions (bar, spheroidals) at different distances from the plane 
(see also Johnson et al. 201 1). This will also enable studies 
of the earliest phases of the Galactic chemical enrichment, 
properly sampling the metal-poor tail of the bulge metallic - 
ity distribution function. It is conjectured that the first stars 
formed are either found in the bulge (Tumlinson 2010) or 
the halo (Sect. 2.3). Since metal-poor stars are very rare in 
the bulge this dissection of the growth and evolution of the 
bulge requires surveys of large numbers of stars (Kunder et 
al. 2012). 

2.3 Galactic halo 

The Galactic halo contains the oldest and most metal-poor 
stars currently known in our Galaxy. Their metallicity distri- 
bution and the abundances of the chemical elements (chem- 
ical abundances) accessible to high-resolution surveys con- 
tain unique imprints of the physical processes that took 
place from the formation of the first stars to the build-up 
of our Galaxy (Karlsson et al. 201 1). In particular, stars be- 
low [Fe/H] ~ — 2.5dex can provide indirect information 
on the star formation processes in metal-poor environments, 
the properties of the first generation of stars (e.g., their ini- 
tial mass function or rotation characteristics), and probe the 
presence of a critical metallicity for star formation (Caffau 
et al. 201 la; Klessen et al. 2012), as well as any dependence 
on environment (as may be the case for different halo build- 
ing blocks). 

However, only three stars at [Fe/H] ^ — 5.0dex are 
known to date. The samples available today are also mostly 
dominated by stars associated with the inner halo, since 



only a small fraction of the samples reach Galactocentric 
distances larger than 15 kpc, where the outer halo popula- 
tion dominates (Carollo et al. 2007). Extrapolating the halo 
metallicity distribution of Schorck et al. (2009), one can ex- 
pect to unearth at least 10 000 objects below [Fe/H]=— 2 dex 
under a typical survey strategy (at some 10 5 targets), still 
yielding several hundreds objects below —4 dex and a few 
tens below —5 dex. 

We also note that Nissen & Schuster (2010) claimed the 
distinction of two groups of stars in the [a/Fe] versus [Fe/H] 
abundance ratio space. These groups can be kinematically 
associated with the inner and outer halo, respectively. Future 
high-resolution spectroscopy halo surveys will extend such 
studies to farther distances, using distances and proper mo- 
tions obtained by Gaia, and applying it to a higher dimen- 
sion abundance ratio space - provided a high level of preci- 
sion in the abundance measurements can be achieved. Like- 
wise, chemical tagging will aid to identify substructures 
(such as hundreds of streams of disrupted satellites) associ- 
ated with past merger and accretion events, as well as char- 
acterise the properties of the progenitor "building blocks", 
such as their chemical enrichment history, and the main nu- 
cleosynthetic channels that led to their present abundance 
patterns (Tolstoy et al. 2009; Koch 2009). 

2.4 Nucleosynthesis 

Detailed chemical abundance distributions will not only 
serve to disentangle the formation of the Galactic compo- 
nents, but can also provide invaluable insight into the origin 
of the elements themselves. Here we briefly touch upon two 
examples that reflect the relevance of high-resolution mea- 
surements of light and heavy elements for the study of nu- 
cleosynthetic processes; their actual measurability will be 
discussed in Sect. 5.2.2. 

According to the standard big bang nucleosynthesis 
model, lithium is produced in primordial nucleosynthe- 
sis. Metal-poor ([M/H]< —1.4) dwarfs, formed from pri- 
mordial material enriched by only a few supernovae, have 
mainly a constant Li abundance, the so called Spite plateau 
dSpite & Spitelll982l) . Recent studies have shown a melt- 
down in this plateau (Sbordone et al. 2010) that is not 
fully understood. As stars evolve to the giant branch, their 
convective zone reaches higher layers in the stellar atmo- 
sphere, and the fragile Li is destroyed by the exposure 
to higher temperatures. Among the targets we plan to ob- 
serve there will be such metal-poor stars, that will allow 
us to enlarge the statistics on the stars on the Spite plateau 



and be yond. It has been recently shown dMucciarelli et al. 



2012albh that old, red giants below the RGB bump dis- 
play a constant lithium abundance that mirrors the Spite 
plateau at a lower level. This observation supports theoreti- 
cal computations that predict that those stars have diluted 
their lithium content, but preserve memory of what their 
original lithium content was. If confirmed by further, ob- 
servations this would open up the possibility of studying 
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the original lithium content in external galaxies in the Lo- 
cal Group, thanks to the higher luminosity of giants with 
respect to dwarfs. 

In spite of constituting only small fraction of the bary- 
onic mass of the MW, n-capture elements leave an impor- 
tant mark on stellar spectra over a large range in metallic - 
ity. The actual formation processes of these heavy elements 
(e.g. the main r-process) remain to date one of the open 
questions in astrophysics. Detailed abundance information 
can allow us to understand the process and the conditions 
under which it works, the sites of the production of these el- 
ements, and their production during different epochs of the 
Galaxy. For a complete review on this subject see Sneden et 
al. (2008). 



2.5 Implications for velocity and abundance 
requirements 

The ambitious goal of the 4MOST survey is to sample 
all major nucleosynthetic channels namely, light elements, 
odd-Z elements (Na, Al), a elements, Fe-peak elements, 
and heavy and light neutron-capture elements such as slow 
neutron-capture (s-) process elements, and rapid neutron- 
capture (r-) process elements. This leads to the requirement 
that we sample at least 15 different elements, as specified in 
Table[U This chemical study of a large sample of Galactic 
stars permits one to better understand our Galaxy. A precise 
knowledge on Fe and a-elements with kinematics would al- 
low one to find and characterise different stellar populations 
in the MW. A chemical map on heavy elements can help 
us understand their actual production processes and define 
the sites of their formation. The precision in chemical abun- 
dances is driven by the need to distinguish different stellar 
populatio ns; for instance, the two halo populations advo- 
cated by Nissen & Schuster (2010) are separated by 0.1- 
0.2 dex in [a/Fe]. 

According to the end-of-mission performances assessed 
at the time of the Gaia Mission Critical Design Review 
(April 201 1 Et), a V=15 mag G2V star will have an error on 
the parallax of 24 /ias. Neglecting interstellar reddening and 
assuming an absolute magnitude of My=4.4 mag for such a 
star, the distance is of the order of 1.3 kpc. This means that 
for a proper motion of 1 mas yr _1 the error on the trans- 
verse velocit}0 is ~ 0.3km s" 1 , for 5 mas yr _1 it is of ~ 
l.lkms -1 . The radial velocities derived from a survey built 
as follow-up of Gaia on the ground should have a compara- 
ble precision. 



3 Targets 

Realistic estimates of spectrograph efficiency, coupled with 
a 4m class telescope lead to a typical limiting magnitude 
of r < 16 mag. The parallaxes from Gaia in this regime will 
permit one to derive distances to within 10% accuracy out 
to ~40 kpc for red giants. The distance scale is, naturally, 
limited by the spectral type and, at V=16mag, F-G dwarf 
distances will be determined at the 10% level out to max- 
imally a few kpc (see also Fig.l in de Jong 2011), which 
bears relevance to contriving the science drivers of follow- 
up surveys (Sect. 2). We concentrate on the case of dwarf 
stars because for them age measurements from isochrones 
will be possible. 

Our science questions above imply that more than 50% 
of the targets observed with the 4MOST high resolution 
mode will be dwarfs with a metallicity in the range —2.0 < 
[Fc/H] < +0.2. Moreover, since the Galactic bulge is 
of crucial interest to the spectroscopic surveys, a fraction 
of the targets will be metal-rich bulge giants. Yet another 
contribution to the samples will come from metal-poor gi- 
ants as found in the Galactic halo. Typical selection crite- 
ria for halo science cases restrict their targets to stars with 
[Fe/H]< -ldex. 

A general concern is that studies of metal-poor giants 
require a predominantly blue wavelength range, while spec- 
troscopy of metal-rich giants benefits from a redder spectral 
coverage. This is due to the fact that for wavelengths below 
about 450 nm, the spectrum of a metal-rich giant is an in- 
extricable series of blends, while for metal-poor giants only 
very few lines are strong enough to be detected for wave- 
lengths red-wards of about 500 nm. 

Chemical abundances for a large number of elements for 
dwarf stars and metal-poor giants can generally be derived 
with high accuracy (see, e.g., Cayrel et al. 2004; Bonifa- 
cio et al. 2009), while this is far more challenging in the 
metal-rich giants of the bulge (Fulbright et al. 2007, Hill 
et al. 2011). Due to the high extinction, typically observ- 
able bulge stars are cool (T e ff < 4200 K). For instance, 



Hill et al. (2011) observed stars in Baade's window with 



the FLAMES/Giraffe instrument at the VLT; their targets 
with an r-magnitude < 16 mag have effective temperatures 
in the range 4300-5300K (see also Johnson et al. 201 1). 

Surveys with 4 m class telescopes are restricted to 
brighter objects, but one will also be able to cover regions 
with higher extinction than the low-extinction windows in 
the Bulge, such as Baade's and Plaut's fields so that typi- 
cal surveys will contain stars with T e fj < 4200 K. For these 
low-temperature stars, the TiO bands in the spectra become 
non-negligible, and spectral synthesis taking into account 
— the molecular bands is mandatory for a d etailed chemical 

3 |http : //www . rssd.esa . int /index .php?project=GATR& pa^g^^ii^^^.^R^<^|^ ottr ^ (2009) list a number of 
« Transverse velocity (V TX proper motion (p) a nd parallax (tt) are j \ { h f f h TiQ b d d f & 

linked through VT= A v ^ jPerrvman & ESAlll99l vol. 1 eqn. 1.2.20), 

a i.ini.n -I-, r i ■ i j • .i fall in the spectral range foreseen for 4MOST (see also Ful- 

Ay = 4.74047 [km yr s J, from which we derive the error expression f o \ 

a yx = A v z*. (0.526 + £) , since ov and a M are correlated through bri g ht et al. 2006). Overall, it is unlikely that the same num- 
a-^ = 0.526<r,r. ber of chemical elements in stars of the bulge (e.g. Bensby 
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Table 1 Number of available lines in the proposed spectral wavelength ranges, separately for dwarfs and giant stars. 



[380; 450 


nm domain 


[580; 680] nm domain 


Parameter 


Element 


# lines 


Element 


# lines 






(dwarf/giant) 


(dwarf/giant) 




?cft 


Fel 


73/70 


Fel 


96/107 








Ha 


1 


log 5 


Fell 


12/7 


Fell 


8/2 


[a/Fe] 


Cal 


8/8 


Cal 


14/14 




Mgl 


2/1 


Mgl 


1/2 




Si I 


3/2 


Sil 


13/12 




Til 


13/19 


Til 


8/27 




Till 


10/26 


Till 


0/1 


[X/Fe] 






Lil 


1 




Nal 


1/0 


Nal 


4/3 




All 


2/1 


All 


1/2 




Scl 


1/0 


Scl 


1/1 




Sen 


3/4 


Sell 


2/6 




VI 


7/7 


VI 


9/27 




VII 


1/1 


VII 


0/0 




Cri 


6/4 


Cri 


1/2 




Cm 


1/1 


Cm 


0/0 




Mm 


19/13 


Mm 


2/3 




Mnll 


0/1 


Mm 


0/3 




Col 


9/11 


Coi 


0/11 




Nil 


3/3 


Nil 


18/19 




Sril 


1/1 


Sril 


0/0 




Yll 


3/1 


Yll 


1/1 




Zrll 


2/2 


Zrll 


0/0 




Ball 


1/1 


Ball 


2/2 




La II 


1/6 


Lai 1 


0/4 




Ndll 


2/5 


Ndll 


1/0 




Eull 


1/1 


Eull 


1/1 




CH 


G-band 








CN 


X-A band 







et al. 2010; Ryde et al. 2010) can be measured, compared 
to, say, subsamples of the Galactic halo (e.g., Fulbright et 
al. 2007; Roederer & Lawler 2012). 

4 Simulated spectra 

We selected stellar parameters that represent typical Galac- 
tic stars and cover the expected parameter space of com- 
mon survey targets. For dwarfs we selected two values of 
effective temperature (6500 K and 5500 K) and two values 
of gravity (logg=3.0 and 4.0 in cgs units). We computed 
synthetic spectra for five metallicities ([M/H] = +0.5, 0.0, 
— 1.0, —2.0, —3.0). For giants we computed one single tem- 
perature (4500 K), two values of gravity (2.0 and 1.0), and 
the same five values of metallicity. The grid contains 30 
models, 20 for dwarfs and 10 for the giants. 

The model grid was computed by using ATLAS9 (Ku- 
rucz 2005) running under Linux (Sbordone 2005). For the 
opacity distribution functions we emplo yed the Solar (PDF) 



and q-enhanced (AODF) versions from lCastelli & Kurucz 



d2003l) . All metal-poor models ([M/H]=-l, -2, -3) have 
been computed with an enhancement in the a-elements 



of 0.4 dex, as is indeed seen in metal poor globular clus- 
ter and halo field stars (e.g., Barbuy 1983; Pritzl et al. 
2005). The micro-turbulence of the models has been fixed 
at 1.0km s _1 , and the mixing-length parameter has been 
adopted as a=1.25. 

We computed the grid of synthetic spectra with the 
SYNTHE code (Kurucz 2005) running under Linux (Sbor- 
done 2005), in the range 370-700 nm. The micro-turbulence 
in the synthesis was set to 1.5kms _1 . No extra rotational 
velocity with respect to the V sim=0.0kms _1 was included 
for all spectral tests. Zero-rotation is a reasonable approx- 
imation for a G-type dwarf, while problems may occur for 
early F-type stars, which can show substantial rotational ve- 
locities (see Gray 2005). In practice, we will derive V sini 
of the stars, either from analysis of the width of the cross- 
correlation peak (Melo et al. 2001), or through Fourier anal- 
ysis of individual line profiles or of the cross-correlation 
function (Diaz et al. 2011). This parameter will then be 
used throughout the chemical analysis. We foresee a stan- 
dard line profile fitting for the abundance analysis. Thus, 
any extra-broadening of the lines due to rotation and ensu- 
ing progressive blending of features will not strongly affect 
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the precision in abundances due to our careful treatment of 
line blends during the fitting process. Only for the cases of 
high rotational velocities in excess of V sim <; 20 kms -1 , 
the analysis will become difficult, but we estimate that the 
fraction of stars showing such properties will be negligible 
in our target samples. 

Likewise, our analysis ignores the effects of micro- 
and macro-turbulence on the uncertainty on the derived 
abundances. These parameters are not physical quantities, 
but fudge-factors introduced in the analysis based on one- 
dimensional, time-independent, stellar model atmospheres. 
The use of hydrodynamical model atmospheres removes the 
need for these parameters (Steffen et al. 2009). By the time 
of the data analysis of the 4MOST data, it can be expected 
that hydrodynamical models will have reached a sufficient 
maturity to be routinely used. Since micro- and macro- 
turbulence depend on the known atmospheric parameters of 
the stars, an adequate calibration can then be reliably de- 
vised (see also Sect. 5.3). 

The original grid has a wavelength sampling corre- 
sponding to a resolving power of 600,000. 

These basic, synthetic spectra were then manipulated to 
provide the input for our tests of observations of resolu- 
tion 15000-25000 as follows. In the first place we de- 
graded the spectra to a resolving power of 250,000 - about 
ten times the value envisioned for a high-resolution spec- 
trograph. Secondly, we normalised the spectral flux to the 
corresponding r-band magnitude as 

F = 3631 * ur- x 10— x (^) 

in [ergs -1 cm -1 A -1 ], where / is the monochromatic flux 
at a reference wavelength of 6166 A, for two reference mag- 
nitudes of r = 14 and 16 mag. Finally, the spectra were re- 
sampled to an equispaced wavelength space with a step size 
in A of 0.5 pm. These spectra have been used to determine 
the spectrograph's design and spectral requirements to reach 
the science goals outlined above. 

The aforementioned synthetic spectra we processed 
through a throughput simulator that produces simulated ob- 
served spectra, taking into account the characteristics of the 
targets (i.e., these very synthetic spectra), a transmission 
model of the Earth's atmosphere, the seeing, the sky back- 
ground, and a simple model of the instrument (see Sartoretti 
et al. 2012 for details). An example of a simulated, noise- 
added spectrum for the case of a G-turn-off star of solar 
metallicity is shown in FigQ] 

The basic characteristics of the spectra computed with 
the simulator that we employ for the following studies are 
reported in Table|2] 

5 High Resolution mode requirements 
5.1 Wavelength range 

For our science case we need an as wide as possible range 
in wavelength. We assume a spectrograph design with two 




580 600 620 640 660 680 

Wavelength [nm] 

Fig. 1 Simulated spectra of the turn-off star 5500/4.0/0.0 
of r=14 mag, observed in New Moon for 6000 s. Top panel 
shows the blue and lower panel the red arm. The S/N ratio 
per pixel at the centre of the blue and red arms is 75 and 
121, respectively. 



Table 2 Spectral characteristics of the simulated spectra. 





Blue arm 


Red arm 


A mi n 


395 nm 


587 nm 




18150 


18635 


Samplingmin 


2pixel/FWHM 


2pixel/FWHM 


Acentre 


425.8 nm 


630 nm 


R-centre 


19563 


20000 


Samplingcentre 


2.51pixel/FWHM 


2.5 pixel/FWHM 


A max 


456.5 nm 


673 nm 


R>max 


20976 


21365 


Sampling max 


4.29pixel/FWHM 


4.22 pixel/FWHM 



arms. If a mosaic of two detectors should be used to cover 
the needed area of the focal plane in each arm, the position 
of the inter-detector gap should also be carefully considered. 

To fix the two wavelength ranges we tried to optimise 
the chemical output, yielding abundances for the largest 
number of elements, but giving more weight to the most 
interesting ones for our science goals. The envisaged wave- 
length ranges are: 

- for the blue arm [392;456.5] nm or [395;459]nm, with 
the provision that the CH G-band at [429;432] nm falls 
completely on one detector; the two possibilities in 
range correspond to the Ca II-K falling inside or outside 
the range, respectively; 

- for the red arm [587;673] nm, with, if two detectors are 
necessary, a gap around 635 nm. 

The above ranges have been selected in order to optimise the 
stellar parameter determination and to derive a maximally 
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comprehensive set of chemical abundance patterns across 
the expected target range. 

An important question concerns the Call-K line 
(393.3663 nm). Call-K is a very strong line and, as such, 
still detectable in the most metal-poor objects. For instance, 
in a star with [M/H] » —5, this line still has an equiva- 
lent width (EW) of about 28 pm, making it too strong to be 
a reliable indicator for the Ca abundance. At a spectral re- 
solving power of R =20 000, this line is thus too strong to 
be useful for the large majority of the stars observed. The 
expected yield of ultra metal poor stars below [M/H]< —5 
is very low (e.g., Christlieb et al. 2004; Frebel et al. 2005; 
Beers & Christlieb 2005; Schorck et al. 2009; Li et al. 2010; 
Caffau et al. 201 la) so that we chose to exclude the K-line 
from the covered wavelength range of the high-resolution 
mode. Ca II-K is a more important feature in low-resolution 
studies (Beers et al. 1985; Beers & Christlieb 2005), where 
lines that strong are the only ones that can still be detected 
at very low metallicities. 

Following common practice, stellar parameters deter- 
minations generally rely on a combination of the indica- 
tors listed below (e.g., Edvardsson et al. 1993; Koch & 
McWilliam 2008; Lee et al. 2008; Sbordone et al. 2010b). 
The cited precisions attainable are what we expect to obtain 
in a survey conducted with such an instrument. 

- T c fj: stems from a combination of photometry, exci- 
tation equilibrium using Fe I lines, and, in F-G dwarf 
stars, from the wings of the Ha line. Temperatures can 
usually be settled to within 150 K or better. 

- log g: employs ionisation equilibrium from Fe I and Fe II 
and, if present, other ionised and neutral species. The 
spectral range indicated above ensures the presence of 
several ionised iron lines. Overall, a precision in the 
gravities within ±0.5 dex is realistic. Ionisation equilib- 
rium of other species can also be used, such as Ti i/Ti II 
or V I/V II. 

- £: micro-turbulent velocities will be determined by 
removing any trend in abundance with line strength 
based on the large number of Fe I lines of various 
strength available in the spectra, to a precision of about 
0.2km s" 1 . 

- [Fe/H]: more than 100 Fe I lines are typically available 
in the envisioned target stars' spectra and our ample 
experience has shown that an accuracy of better than 
0. 1 dex is easily achievable. 

- [a/Fe]: several Ca and Mg lines are measurable in the 
spectra (see TableQ} and precisions of 0. 1 dex are com- 
monly reached. 

- [X/H] for a variety of chemical species X: ideally, 
spectroscopic surveys foresee precisions reaching from 
0.1 dex (e.g., in [Ti/H]) to 0.5 dex (e.g., in [S/H], where 
feasible). 

Given these constraints, we summarise the useful (i.e., 
unblended and/or recoverable) lines to be employed in a 
meaningful, high-resolution abundance analysis in TableQ] 



Here, we list the number of lines that are potentially de- 
tectable in the selected wavelength ranges, separately for 
dwarfs and giants, although we caution that not all lines will 
be detectable in all spectra due to their different sensitivities 
to the stellar parameters. Hence, in a typical survey strategy, 
for those elements, for which only few lines are available, 
one cannot derive a chemical abundance ratio for the com- 
plete set of targeted stars. The spectral content, summarised 
in TableQ] contains diagnostic information for the determi- 
nation of the atmospheric parameters. 

In addition to the aforementioned Fe- and a-element 
transitions, we will rely on the following absorption features 
for the abundance measurements of further, crucial tracers 
of Galactic chemical evolution: the G-band of CH (429- 
432 nm), to derive the abundance of carbon; the CN band at 
414^422 nm will be used as an indicator of the N abundance 
once the C-abundance is known; the [OI] 630 nm feature 
can serve to derive O-abundances in giants; the Nal-lines 
at 589 nm and All lines (at 396. 152 nm) cover representa- 
tive p-capture elements; the Li I doublet at 670.7 nm is the 
only accessible means of determining A(Li). Finally, sev- 
eral transitions are available to derive the abundance ratios 
for some n-capture elements. Sr, Y, and Zr trace the first 
peak of n-capture elements and Ba, La, Ce, and Nd trace 
the second peak. In the solar system abundances, all have 
been produced mainly through the slow n-capture process 
(s-process), although all have a rapid n-capture process (re- 
process) component that can be as large as 42% in the case 
of Nd (Sneden et al. 2008). Finally also Eu, an almost pure 
r-process element, can be observed in the proposed 4MOST 
configuration. 

We do not consider in detail the G-band in our further 
tests, since this feature is clearly measurable also at low res- 
olution, and thus would put no stringent constraint on the 
design of the HR spectrograph. 

We finally mention some interesting features that one 
looses with our trade-off choice of the wavelength ranges. 

- The Mgl-b triplet at about 516nm, which is useful for 
abundance determination in metal-poor stars, [M/H] < 
—2.5, and also for gravity determination for metal-rich 
stars. The large majority of stars that 4MOST will ob- 
serve will be close to solar metallicity, and, at the resolv- 
ing power of 4MOST, the lines of the Mg b triplet are 
strongly saturated, rendering them not useful for abun- 
dance determination. The wings of the saturated Mgb 
triplet lines would be a useful indicator for the stellar 
gravity, but our spectral range already contains a suffi- 
cient number of gravity-sensitive features (TableQ]i and 
the majority of the stars, as stated above, will be selected 
to be turn-off stars, so that their gravity will be known a 
priori. 

- The oxygen triplet at 777 nm is a very useful tool to de- 
rive the oxygen abundance, in dwarf stars at all metallic- 
ities. Oxygen abundance in giants will be derived from 
the [OI] line at 630 nm. For dwarf stars we will not be 
able to measure any O-abundance. While including the 
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oxygen triplet would thus be highly valuable, the range 
around it is very poor in lines so that no other informa- 
tion could be gleaned. 
- The near-infrared Call triplet at about 850 nm is in- 
cluded in the range observed by Gaia. We excluded 
these features from our HR wavelength ranges for the 
same reasons as for excluding the Mg b triplet and oxy- 
gen triplets, i.e., the lines are too strong to be good abun- 
dance indicators in solar metallicity stars at the resolv- 
ing power of 4MOST and the range around them has lit- 
tle other interesting lines. Furthermore, that wavelength 
range is highly contaminated by telluric absorption and 
night-sky emission lines. 



5.2 Spectral resolution and pixel-size 

In the design of a spectrograph, resolution and spectral cov- 
erage are correlated. In order to avoid undersampling one 
needs at least two pixels per resolution element. This im- 
plies that for a detector of given physical dimension of the 
pixel and number of pixels increasing the resolution reduces 
the spectral coverage. While a higher resolution is desirable 
for a chemical abundance analysis, we also prefer a wave- 
length range that is as wide as possible. This should be cou- 
pled to an overall instrumental efficiency that allows for the 
detection of the spectral lines over the spectral noise. 

Another aspect that influences the chemical line-by-line 
analysis is pixel-size (in wavelength units, defined as pixel- 
size= Aj+i — Aj, see Cayrel 1988): a small pixel-size al- 
lows one to better characterise the shape of each line that 
is sampled by a large number of pixels. However, a larger 
number of pixels decreases the S/N ratio per pixel, be- 
cause the same number of photons is spread over more pix- 
els. An unblended line, without hyperfine structure and/or 
isotopic splitting, needs several pixels (from experience at 
least about five, one around the line core and two on each 
wing) for its shape to be adequately described, in a way to 
allow line profile fitting for the abundance analysis. The full 
line width depends on its intensity (see appendix A) so that 
weaker lines are sampled by a smaller number of pixels. 
Thus the smallest acceptable pixel-size is fixed by the weak- 
est lines that we want to measure. For a fixed resolution the 
spectral sampling (i.e., the number of pixels per resolution 
element) changes with wavelength: 
sampling = A/(resolving power)/pixel-size. 

It is self-evident that the instrumental characteristics 
must be feasible from a technical point of view. Therefore, 
we investigate in the following the cases of a good compro- 
mise, at an average resolving power of 20 000 and a min- 
imum requirement of R > 15 000 across the entire spec- 
tral range, and a sampling between 2 and about 4, which is 
sufficient to resolve the lines without dividing the photons 
between too many pixels. 

www.an-journal.org 



5.2.1 Weak versus strong lines 

As a first test, we investigated two cases: the strong Ba II 
line at 455.4nm (6pm < EW < 18pm) and the Ti I line at 
454.8 nm as an example of a weak line (EW=2.3 pm for the 
metal-poor model with T off /log g/[MM] = 5500/4.0/-2.0). 

The Ball line lies, especially in metal-rich stars, in 
a crowded spectral region. Thus this feature allows us to 
probe also the effects of blending. In practice, we fitted, 
minimising the x 2 , the line profile in a synthetic spectrum 
that has been injected with Poisson noise to mimic S/N 
ratios of 50 and 40 per pixel at three different values for re- 
solving powers 15 000, 20000, and 25 000 and a pixel-size 
fixed at 2 pm. The results are reported in Table[3] where we 
list the recovered, mean abundance ratios and their standard 
deviations (a). The latter were obtained by running 1 000 
Monte Carlo trials. In some cases the fit was not success- 
ful, giving an abundance which deviates from the input one 
by several orders of magnitudes, and we reject these events, 
while keeping all others that we define as "good events". As 
expected, as the resolution decreases, the numbers of times 
the fit fails increases, as does the scatter a. 

In the next step, we investigated the effect that the com- 
bination of resolving power with a given pixel-size has on 
the abundance determination. For each of the above resolv- 
ing powers, we considered different samplings of 2.0, 2.5 
and 3.0 pixels per resolution element. The resulting pixel- 
size was computed at a reference wavelength of 395 nm. 
In general, very strong lines (EW > 10 pm) are easily 
detectable above spectral noise even at faint magnitudes. 
However, such features are less suitable for a reliable abun- 
dance determination, because they are difficult to model due 
to uncertain damping constants and usually stronger depar- 
tures from Local Thermodynamic Equilibrium (LTE), and 
the growing importance of 3D effects (e.g., Caffau et al. 
2011b). Furthermore, they are significantly less sensitive 
to changes in abundance compared to weaker, unsaturated 
lines. The EWs quoted in Table|4]refer to the suitably strong 
Ba II transition in the metal-poor model and the weaker Ti I 
line in the same model spectrum. 

As these tests indicate, a strong line is much less af- 
fected by a change in resolving power than a weak line. In 
contrast, changes in pixel-size at a given S/N ratio (stated 
per A) have no significant effect on the abundance deter- 
mination. A consequence of fixing the S/N ratio per A, 
a decrease in pixel-size leads to a better sampling of the 
line profile, with a higher number of pixels that resolve the 
line, but at the expense of a lower S/N ratio in each pixel. 
In conclusion, we suggest as a compromise in resolving 
power and pixel sampling one of the following combina- 
tions: 20 000/2.5, 25 000/2.0-2.5, or 15 000/2.5-3.0. 

5.2.2 Particular lines 

In the following we briefly reflect upon a few elements 
of particular importance (Sect. 2.4), for which abundances 
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Table 3 Accuracy of the 455.4 nm Ball line as a function of resolution. The pixel-size is 0.002 nm. The S/N is given 
per pixel. 



Model (T Pff /log ff/[Fe/H] ) 


Resolving Power S/N 


Good events 


[Ba/H] 




EW = 


17.5 pm 






4500/2.0/- 1.0 


15000 


50 


976/1000 


-1.021 ±0.065 


4500/2.0/- 1.0 


20000 


50 


1000/1000 


-1.027 ±0.050 


4500/2.0/- 1.0 


25 000 


50 


1000/1000 


-1.029 ±0.047 


4500/2.0/- 1.0 


15000 


40 


972/1000 


-1.017 ±0.079 


4500/2.0/- 1.0 


20000 


40 


1000/1000 


-1.023 ±0.061 


4500/2.0/- 1.0 


25 000 


40 


1000/1000 


-1.026 ±0.056 




EW = 


17.0 pm 






5500/4.0/+0.0 


15000 


50 


996/1000 


-0.013 ±0.064 


5500/4.0/+0.0 


20000 


50 


999/1000 


-0.012 ±0.048 


5500/4.0/+0.0 


25 000 


50 


1000/1000 


-0.019 ±0.039 




EW = 


1 1.6 pm 






5500/4.0/- 1.0 


15000 


50 


974/1000 


-0.995 ±0.071 


5500/4.0/- 1.0 


20000 


50 


997/1000 


-1.003 ±0.051 


5500/4.0/- 1.0 


25 000 


50 


1000/1000 


-0.997 ±0.044 



Table 4 Accuracy of the 455.4 nm Ba II line as a function of resolving power and pixel size. These tests are for a strong 
line and the 454.8 nm Ti I line as a test case for a weak line. The S/N is 160 per A each. 



Line 


Model 




Pixel 


Resolving Power 


S/N 


Sampling 


Good 


[X/H] 


[nm] 


(T cff /log ff /[Fe/H]) 


Size [nm] 




per pixel 


@395nm 


Events 




EW = 6.6 pm 


Ball 455.4 


5500/4.0/- 


2.0 


0.0078 


25000 


45 


2.0 


977/1000 


-1.965 ±0.098 


Ba II 455.4 


5500/4.0/- 


2.0 


0.00624 


25 000 


40 


2.5 


976/1000 


-1.969 ±0.093 


Ball 455.4 


5500/4.0/- 


2.0 


0.0052 


25 000 


36 


3.0 


985/1000 


-1.973 ±0.091 


Ball 455.4 


5500/4.0/- 


2.0 


0.00975 


20000 


50 


2.0 


979/1000 


-1.959 ±0.113 


Ba II 455.4 


5500/4.0/- 


2.0 


0.0078 


20000 


45 


2.5 


971/1000 


-1.964 ±0.115 


Ba II 455.4 


5500/4.0/- 


2.0 


0.0065 


20000 


41 


3.0 


973/1000 


-1.970 ±0.111 


Ball 455.4 


5500/4.0/- 


2.0 


0.0130 


15 000 


58 


2.0 


929/1000 


-1.950 ±0.157 


Ba II 455.4 


5500/4.0/- 


2.0 


0.0104 


15 000 


52 


2.5 


938/1000 


-1.950 ±0.162 


Ball 455.4 


5500/4.0/- 


2.0 


0.0087 


15 000 


47 


3.0 


903/1000 


-1.953 ±0.162 


EW = 2.0 pm 


Ti I 454.8 


5500/4.0/- 


2.0 


0.0078 


25 000 


45 


2.0 


814/1000 


-2.002 ±0.137 


Ti I 454.8 


5500/4.0/- 


2.0 


0.00624 


25 000 


40 


2.5 


763/1000 


-2.000 ±0.150 


Ti I 454.8 


5500/4.0/- 


2.0 


0.0052 


25 000 


36 


3.0 


836/1000 


-2.011 ±0.142 


Ti I 454.8 


5500/4.0/- 


2.0 


0.00975 


20000 


50 


2.0 


741/1000 


-2.154 ±0.237 


Ti 1 454.8 


5500/4.0/- 


2.0 


0.0078 


20000 


45 


2.5 


764/1000 


-2.008 ± 0.180 


Ti I 454.8 


5500/4.0/- 


2.0 


0.0065 


20000 


41 


3.0 


640/1000 


-2.027 ±0.202 


Ti I 454.8 


5500/4.0/- 


2.0 


0.0130 


15 000 


58 


2.0 


645/1000 


-1.973 ±0.301 


Ti I 454.8 


5500/4.0/- 


2.0 


0.0104 


15 000 


52 


2.5 


651/1000 


-1.973 ±0.258 


Ti I 454.8 


5500/4.0/- 


2.0 


0.0087 


15 000 


47 


3.0 


597/1000 


-1.979 ±0.315 



need to be derived from essentially one line and which re- 
quire special care. 

Oxygen is a fundamental tracer of fast chemical en- 
richment that can efficiently be derived in giants, but poses 
several problems in dwarf spectra. The forbidden [OI] line 
at 630 nm (see Table[5]), available in our wavelength range 
(and a typical EW of about 0.35 pm in the solar spectrum; 
Caffau et al. 2008) is too weak in dwarf stars. The line 
becomes stronger in giants and can be measured down to 
[Fe/H]< —3 (e.g., Cayrel et al. 2004). Some telluric absorp- 
tions affect the region, and, for stars with low radial veloc- 
ity, the atmospheric emission can also affect the line. Even 



if the line is problematic, it is very important to derive oxy- 
gen abundance, at least for a fraction of stars. As our tests 
show, a resolving power of at least 20 000 is necessary to de- 
rive the oxygen abundance from this [OI] line, see Table[5] 
In fact, our simulations indicate that a resolving power of 
15 000 yields only 60% "good events" (for the case of the 
model at [M/H]=— 1.0), compared to a success rate of 90% 
at a resolving power of 25 000. 

The lithium abundance is mostly derived from the 
670.7 nm doublet. Our experience has shown that this fea- 
ture can be detected when the spectral resolving power is 
20000, provided a S/N ratio of about 80 per pixel (see 
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Table 5 The same as Tableg] but for the [OI] line at 630.0 nm. 



Line 


Model 


Pixel 


Resolving Power 


S/N 


Sampling 


Good 


[O/H] 


[nm] 


(T eff /logg/[Fe/H]) 


Size [nm] 






per pixel 


@395 nm 


Events 












EW = 


4.8 pm 










[OI] 


630.0 


4500/2.0/0.0 


0.0078 


25 000 




45 


2.0 


919/1000 


0.013 ±0.104 


[OI] 


630.0 


4500/2.0/0.0 


0.00624 


25 000 




40 


2.5 


920/1000 


0.019 ± 0.108 


[OI] 


630.0 


4500/2.0/0.0 


0.0052 


25 000 




36 


3.0 


942/1000 


0.016 ± 0.108 


[OI] 


630.0 


4500/2.0/0.0 


0.00975 


20000 




50 


2.0 


905/1000 


0.015 ±0.144 


[OI] 


630.0 


4500/2.0/0.0 


0.0078 


20000 




45 


2.5 


904/1000 


0.020 ±0.133 


[OI] 


630.0 


4500/2.0/0.0 


0.0065 


20000 




41 


3.0 


920/1000 


0.026 ± 0.133 


[OI] 


630.0 


4500/2.0/0.0 


0.0130 


15 000 




58 


2.0 


893/1000 


0.036 ±0.198 


[OI] 


630.0 


4500/2.0/0.0 


0.0104 


15 000 




52 


2.5 


892/1000 


0.042 ±0.182 


[OI] 


630.0 


4500/2.0/0.0 


0.0087 


15 000 




47 


3.0 


882/1000 


0.038 ± 0.201 










EW = 


0.8 pm 










[OI] 


630.0 


4500/2.0/- 1.0 


0.0078 


Zj uuu 




45 


2.0 


792/1000 


-0.974 ±0.185 


[OI] 


630.0 


a ^nn/o n/ i n 
4jUU/z.U/-1.U 


U.UUDZ4 


25 000 




-in 


Z.J 


Q/ic/i nnn 
o4j/1UUU 


ft 07Q J-fl 1 QQ 


[OI] 


630.0 


4500/2.0/- 1.0 


0.0052 


25 000 




36 


3.0 


838/1000 


-0.964 ±0.180 


[OI] 


630.0 


4500/2.0/- 1.0 


0.00975 


20000 




50 


2.0 


749/1000 


-0.941 ± 0.220 


[OI] 


630.0 


4500/2.0/- 1.0 


0.0078 


20000 




45 


2.5 


695/1000 


-0.967 ±0.231 


[OI] 


630.0 


4500/2.0/- 1.0 


0.0065 


20000 




41 


3.0 


759/1000 


-0.940 ±0.215 


[OI] 


630.0 


4500/2.0/- 1.0 


0.0130 


15 000 




58 


2.0 


609/1000 


-0.790 ± 0.329 


[OI] 


630.0 


4500/2.0/- 1.0 


0.0104 


15 000 




52 


2.5 


623/1000 


-0.778 ± 0.349 


[OI] 


630.0 


4500/2.0/- 1.0 


0.0087 


15 000 




47 


3.0 


590/1000 


-0.772 ± 0.339 



e.g. Monaco et al. 2010, 2012). With the spectrograph de- 
sign we are testing in this work, this S/N ratio is achieved 
only for the brightest objects. Nevertheless, this will be suf- 
ficient to derive Li abundance for a large sample of Galactic 
stars to allow us to improve our understanding on the history 
of lithium in the MW. 

Strontium is the first of the weak s-process neutron- 
capture elements at Solar metallicity (Heil et al. 2009, Pig- 
natari et al. 2008, and Arlandini et al. 1999), which we can 
detect. Due to the strong resonance lines (in particular the 
407.7 nm line) the detection of this element sets no con- 
straint on the required resolution. It remains measurable in 
both dwarfs and giants at all metallicities. However, in the 
stars [Fe/H] > —2.0 the measurement of this line poses a 
problem in an automatic fitting routine owing to the strong 
line blends in both Sr wings, and the fact that the line satu- 
rates at even higher metallicities. 

As an example for a rapid n-capture element, we consid- 
ered Eu and ran our tests on the Eu II line at 664.5 nm (see 
Tables[T]and|6]l. Europium is a second peak neutron-capture 
element. It is mainly synthesised by r-process in massive 
stars (not yet fully understood - for details see Francois et 
al. 2007). As our tests indicate, this procedure yields reli- 
able results both for dwarf and giant stars so that no strin- 
gent limitations in resolution is required in this case. 

5.3 Precision in chemical abundances 

We would like to stress that the precision we are here con- 
sidering is only referring to the random errors based on the 
noise characteristics of the spectra. No systematic error re- 
lated to uncertainties on stellar parameters (T e ff, logg, £, 



[M/H]) and atomic data (^/-values) will be taken into ac- 
count, neither other noise sources. 

With millions of individual stars, the amount of spec- 
tral data that present-day surveys wi ll observe is too large 



for m anual abundance analysis (e.g. iBonifacio et alj 2009 



Cavrel et al.ll2004l) . Thus, an automatic approach for deriv- 



ing the stellar parameters and the chemical compositions of 
the observed stars is necessary. Only a few species (e.g., Li, 
S, Sr, and Th) require an ad hoc analysis approach. For in- 
stance, Li requires a NLTE treatment, Sr is heavily blended, 
and the Th line at 408.6 nm lies in the red wing of a stronger 
blend that has to be taken into account properly. 

In the present analysis, we used the code MyGIsFOS 
(Sbordone et al. 2010a) to derive the chemical composition 
of our test cases. MyGIsFOS performs a fully automated 
abundance analysis. From a number of predefined regions, 
MyGIsFOS performs a "pseudo-normalisation" of both the 
input spectrum to be analysed and the grid of synthetic spec- 
tra. The chemical abundance analysis is then performed by 
line profile fitting of a sample of selected lines. In Fig. [2] 
an example of the analysis of MyGIsFOS is shown. One of 
the advantages of the line-profile-fitting technique is that, 
in principle, one is free from uncertainties linked to blends. 
This is true in the present analysis since the atomic data of 
the simulations are exactly the same as the ones of the syn- 
thetic grid. In reality, the blends introduce an uncertainty 
due to poor knowledge of the atomic data of the blending 
lines and the potential presence of unknown blends. An as- 
sessment of this kind of uncertainty is beyond the scope of 
the present analysis. 

With MyGIsFOS we analysed the spectra obtained from 
the 4MOST simulator, computed with the characteristics 
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Table 6 The same as Table|5] but for the Eu II line at 664.5 nm. 



Line 


Model 


Pixel 


Resolving Power 


S/N 


Sampling 


Good 


[Eu/H] 


[nm] 


(T cff /log ff /[Fe/H]) 


Size [nm] 




per pixel 


@395nm 


Events 




EW = 5.0 pm 


Eu II 664.5 


4500/2.0/0.0 


0.0078 


25000 


45 


2.0 


901/1000 


0.000 ±0.101 


Eu II 664.5 


4500/2.0/0.0 


0.00624 


25 000 


40 


2.5 


896/1000 


0.001 ±0.101 


Eu II 664.5 


4500/2.0/0.0 


0.0052 


25 000 


36 


3.0 


913/1000 


0.002 ±0.108 


Eu II 664.5 


4500/2.0/0.0 


0.00975 


20000 


50 


2.0 


915/1000 


0.001 ±0.116 


Eu II 664.5 


4500/2.0/0.0 


0.0078 


20000 


45 


2.5 


927/1000 


0.001 ±0.118 


Eu II 664.5 


4500/2.0/0.0 


0.0065 


20000 


41 


3.0 


906/1000 


-0.001 ±0.111 


Eu II 664.5 


4500/2.0/0.0 


0.0130 


15 000 


58 


2.0 


836/1000 


0.028 ±0.136 


Eu II 664.5 


4500/2.0/0.0 


0.0104 


15 000 


52 


2.5 


892/1000 


0.014 ±0.136 


Eu II 664.5 


4500/2.0/0.0 


0.0087 


15 000 


47 


3.0 


882/1000 


0.014 ±0.139 


EW = 10.0 pm 


Eu II 412.9 


4500/2.0/-2.0 


0.0078 


25 000 


45 


2.0 


952/1000 


-1.941 ±0.394 


Eu II 412.9 


4500/2.0/-2.0 


0.00624 


25 000 


40 


2.5 


945/1000 


-1.951 ±0.337 


Eu II 412.9 


4500/2.0/-2.0 


0.0052 


25 000 


36 


3.0 


969/1000 


-1.961 ±0.302 


Eu II 412.9 


4500/2.0/-2.0 


0.00975 


20000 


50 


2.0 


900/1000 


-1.954 ±0.273 


Eu II 412.9 


4500/2.0/-2.0 


0.0078 


20000 


45 


2.5 


953/1000 


-1.956 ±0.311 


Eu II 412.9 


4500/2.0/-2.0 


0.0065 


20000 


41 


3.0 


958/1000 


-1.942 ±0.365 


Eu II 412.9 


4500/2.0/-2.0 


0.0130 


15000 


58 


2.0 


823/1000 


-1.893 ±0.554 


EuH 412.9 


4500/2.0/-2.0 


0.0104 


15000 


52 


2.5 


931/1000 


-1.936 ±0.393 


Eu II 412.9 


4500/2.0/-2.0 


0.0087 


15 000 


47 


3.0 


933/1000 


-1.936 ±0.400 



given in Table[2] (the overall transmission of the instru- 
ment plus telescope is shown in Figj3) to derive the de- 
tailed chemical abundances, while temperature, gravity, and 
micro-turbulence are kept fixed at the value of the input sim- 
ulated synthetic spectrum. 




s _ 

o 
o 

p - 

H 1 1 1 1 1 1 H 

350 400 450 500 550 600 650 700 
wavelength (nm) 

Fig. 3 The transmission of the instrument plus the tele- 
scope. 

As an example of our procedure, we report in Tables|7] 
and[8]the chemical abundances we derive for a Solar metal- 
licity turn-off star (5500/4.0/0.0) and a metal-poor giant 



(4500/2. 0/-3.0), both at an adopted magnitude of r = 
14 mag and simulated to be observed for 3 200 s, which is 
a representative exposure time in 4 m class telescope survey 
strategies. 

The turn-off star of Solar metallicity is a typical star of 
the Galactic disc. The line-to-line scatters^, a, in Table|7] 
are mostly within 0. 1 dex for all elements, except for Zr, for 
which the abundance determination would greatly benefit 
from an increase in exposure time. 

For the Galactic halo star (Tabled the considerations 
are necessarily different. The average metallicity of inner 
halo stars is —1.6 dex (Carollo et al. 2007). When observing 
giants, one can observe stars further away from the Sun and 
thereby probe the halo star composition. One of the goals 
in the target selection of halo stars is to detect and charac- 
terise the most metal-poor stars in the Galactic halo. The 
lower the content of metals in a star the weaker the metallic 
lines appear in the spectrum. Thus we decided to test as an 
extreme case the most metal-poor giant, for which we com- 
puted a synthetic spectrum with [M/H] = —3.0. The results 
will then represent a lower limit to the precision compared 
to what can be achieved for stars 10-100 times more metal- 
rich. 

The line-to-line scatter listed in Tables|7]and|8j both for 
the solar turn-off and for the metal -poor giant stars, is within 
0.1 dex for the large majority of the elements. It becomes 
significantly larger than 0. 1 dex for Zr in the turn-off spec- 
trum, where there are only two lines of Zr II. The same holds 
for Si, Sc, and La for the giant. For some elements a sys- 
tematic effect is visible, in that the average abundance de- 

5 That is, the standard deviation of the abundance measurements from 
all lines of one element 
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Fig. 2 A portion of the spectrum of a turn-off star with Solar-scaled abundances. Gray boxes show the pseudo-continuum 
window, while coloured boxes mark windows around several absorption features: Fe I (blue), a-elements (green), and 
heavy elements (red). Here, the red boxes show the fits of the Ball line at 455.4 nm and of the Cm feature at 455.8 nm. 
The gray shaded line is rejected because the \ 2 is bellow the fixed threshold. Magenta horizontal lines: solid lines represent 
the determined continuum, dashed lines represent one a of the S/N ratio, dotted lines represent three a of the S/N ratio. 
Magenta vertical lines show the shift performed by the fitting procedure. 



rived from the analysis differs from the known input value. 
This effect is more evident for elements whose abundance 
is based on a single line. 

A similar situation holds for the faintest targets at r = 
16 mag, assuming an exposure time of 6800 s. This results 
in a small increase of the uncertainties. We analysed also 
a metal-poor turn-off star (5500/4.0/— 2.0), that shows in 
its spectrum lines weaker than the same lines in the solar- 
metallicity spectrum. As expected, the uncertainties we de- 
rive (see Table|7]i are larger than in the solar-metallicity case 
for some elements 

5.4 Monte Carlo simulation on the complete 
wavelength range 

In the previous sections, we analysed single simulated spec- 
tra of individual targets, computed by the 4MOST simula- 
tor. A single case of noise-injection can be useful to ob- 
tain a first, overall characterisation of the performance of a 
spectrograph under a fixed design, e.g. the minimum EW 
measurable in a line as a function of the S/N ratio, but it is 
generally not sufficient to obtain a careful understanding on 
the final precision achievable, especially for those elements 
where abundances are based on a few lines only (Table[T|i. 

Therefore, with the 4MOST simulator we computed a 
series of 1 000 simulations of the same spectrum, based on 
the values reported in Table|2] These were analysed with 
a Monte Carlo technique and employing the full range in 



wavelength to get an estimate of the overall performance of 
a high resolution spectroscopic survey. In real astronomical 
observations it is customary to estimate the precision on the 
abundance of a given element from the line-to-line scatter. 
The Monte Carlo simulation provides the mean line-to-line 
scatter and its dispersion, over the sample, for all elements 
whose abundance can be derived from several lines. For all 
elements the dispersion on the abundance, over the sam- 
ple is also available. We consider the former, when about 
ten lines of the element are measurable, as an estimate of 
the achievable precision, resorting to the latter for elements 
measurable from a few lines. The two estimates are clearly 
connected, the latter being roughly the former divided by 
the square root of the number of lines. 

The results are shown in Table|9] where we report the 
results for a solar metallicity turn-off star, and in TablefTol 
where we show the recovered abundances for a metal-poor 
turn-off star. In the latter case, Cr, Eu, Nd, and ionised Si, 
present in the solar metallicity case, cannot be analysed, 
since the lines become too weak. Vice versa, Sr is not mea- 
surable in an automated way in the solar metallicity star, 
due to the excessive blending, while it can be measured au- 
tomatically in the low metallicity case. Generally, the aver- 
age line-to-line scatter is about a factor of two larger for the 
metal-poor case, but there are exceptions, such as Ba, whose 
line-to-line scatter is smaller in the metal-poor model, or 
Fe II, for which both models show the same value. 
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Table 7 Chemical abundance results for a solar metallic- 
ity turn-off star (5500 K/4.0/0.0), with r = 14 mag. a de- 
notes the line-to-line scatter. 



Element 


N lines 


[X/H] 


a 


Nal 


2 


-0.10 


0.02 


Mgl 


2 


-0.01 


0.06 


All 


2 


-0.13 


0.03 


Si I 


11 


-0.06 


0.04 


Cal 


16 


-0.03 


0.04 


Sc I/Sc II 


2/4 


+0.07/-0.10 


0.11/0.07 


Til/Till 


18/6 


-0.04/-0.04 


0.07/0.06 


VIA' II 


13/1 


-0.04/-0.25 


0.07 


Cri/Crll 


3/1 


+0.05/-0.011 


0.03 


Mm 


13 


-0.11 


0.03 


Fe I/Fe I 


116/11 


-0.07/-0.12 


0.06/0.09 


Col 


6 


-0.07 


0.07 


Nil 


16 


-0.05 


0.05 


Yll 


2 


-0.10 


0.02 


Zrll 


2 


+0.14 


0.35 


Ban 


2 


-0.13 


0.02 


La II 


1 


-0.08 




Ndn 


2 


-0.08 


0.07 


Euli 


1 


-0.14 





Table 8 The same as Table|7] but for a giant star 
(4500K/2.0/-3.0), with r = 14 mag. 



Element 


N lines 


[X/H] 




a 


Nal 


2 


-3.13 




0.03 


Mgl 


1 


-2.55 






All 


1 


-2.99 






Si I 


2 


-2.42 




0.20 


Cal 


13 


-2.60 




0.04 


Sell 


4 


-3.11 




0.20 


Ti I/Ti II 


15/19 


-2.58/- 


-2.57 


0.08/0.07 


VI 


3 


-2.96 




0.02 


Crl 


3 


-2.97 




0.05 


Mn I/Mn II 


7/1 


-3.00/- 


-2.93 


0.05 


Fe I/Fe II 


63/4 


-2.99/- 


-2.84 


0.07/0.11 


Col 


4 


-3.06 




0.07 


Nil 


4 


-2.88 




0.10 


Sril 


1 


-3.08 






Yll 


1 


-3.08 






Zrll 


2 


-3.01 






Ball 


2 


-3.02 




0.07 


La II 


5 


-2.92 




0.23 


Ndll 


1 


-2.93 






Euli 


1 


-3.11 







As a result, in the case of the Solar turn-off star, we re- 
cover the abundance ratios of the majority of the elements 
systematically lower by typically a few hundredths of dex 
compared to the input spectrum of known parameters and 
input abundances. We also ran the Monte Carlo tests for 
the metal-poor giant with two different values of exposure 
time, 3 200 and 6 000 s (see TablefTTTi, and for the solar- 
metallicity giant (see Tablell2li. As expected, by increasing 
the integration time, the scatter of the average value de- 



Table 9 Abundances results and uncertainties derived 
from 1 000 Monte-Carlo simulations for the (5500/4.0/0.0) 
model with r = 14 mag. In the first column the elements 
analysed are listed. The following three columns refers to 
1 000 Monte Carlo simulation events with an exposure time 
of 3 200 s, and are the minimum-maximum number of lines 
used for the abundance determination, the average abun- 
dance over the 1 000 Monte Carlo events, with its scatter, 
and the average line-to-line scatter over the 1 000 events 
with related standard deviation. 



Element 


N lines 


(fX /HI) 




(7 




Na I 


2 


— 0.04:2 ± 0.058 


0. 


.026 ± 


.020 


Mgl 


1-2 


r\ /"i i~i / * i n nor 

0.026 ± 0.035 


0. 


.042 ± 0. 


.031 


Al I 


1-2 


0.038 ± 0.055 


0. 


.050 ± 


.038 


c : t 
Si I 


6-14 


—0.038 ± 0.027 


0. 


.054 ± 0. 


.014 


Sill 


1-2 


0.014 ± 0.074 


0. 


.083 ± 0. 


.064 


Cal 


8-16 


— 0.017 ± 0.022 


0. 


.045 ± 


.009 


Sen 


2-5 


-0.010 ±0.055 


0. 


.080 ± 


.033 


Til 


8-18 


-0.036 ±0.026 


0. 


.060 ± 


.014 


Till 


2-7 


0.018 ±0.041 


0. 


.057 ±0 


.028 


VI 


7-14 


-0.043 ± 0.021 


0. 


.062 ± 


.014 


Crl 


2-A 


-0.034 ± 0.048 


0. 


.042 ± 


.026 


Cm 


0-1 


0.028 ±0.108 








Mm 


3-12 


-0.064 ±0.039 


0. 


.061 ±0 


.021 


Fel 


98-128 


-0.056 ±0.009 


0. 


.070 ± 


.005 


Fe II 


7-16 


-0.050 ± 0.036 


0. 


.101 ±0 


.020 


Col 


2-6 


-0.089 ± 0.037 


0. 


.039 ± 


.028 


Nil 


8-20 


-0.056 ± 0.024 


0. 


.051 ±0 


.011 


Yll 


1-2 


-0.031 ± 0.068 


0. 


.048 ± 


.046 


Zrll 


1-2 


-0.013 ± 0.067 


0. 


.071 ±0 


.065 


Ball 


0-2 


0.014 ±0.085 





.086 ± 


.068 


La II 


0-1 


-0.082 ±0.111 








Ndll 


1-2 


-0.036 ± 0.091 





.080 ± 


.078 


Euli 


0-1 


-0.023 ± 0.048 









creases, as well as the average line-to-line scatter. Twenty 
elements can be detected in the halo giant, and for five ele- 
ments we can derive the ionisation equilibrium. 

5.4.1 Precision on abundances as a function of 
signal-to-noise ratio 

In the next step, we tested the precision we can achieve 
in the chemical abundance determination as a function of 
the S/N ratio. We concentrated on the solar dwarf case 
(5500/4.0/0.0, r = 14 mag), and run 1000 Monte Carlo 
simulations with six values of exposure time (1 200, 2400, 
3 200, 4 800, 6 000, and 7 200 s). The S/N ratios per pixel, 
for three different wavelengths on either CCD, in the case of 
exposure time of 3 200 s for r = 14, are given in Table[T3l 
The results are summarised in Fig.|4]and|5]for the cases of 
Fe, Ca, Ni (Fig. 2), and Zr, Ba, Eu (Fig. 3), respectively. 

For the elements in Fig. [4] the abundance is based on 
several lines, so that the line-to-line scatter has a meaning. 
For the elements in Fig. [5] only very few lines are used for 
the abundance determination, so that the scatter on the aver- 
age value over the 1 000 Monte Carlo events is a more robust 
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Table 11 Abundances results and uncertainties derived from 1 000 Monte-Carlo simulations for the (4500/2.0/— 3.0) 
model and r = 14 mag. Columns are as in Table|9] but there are two Monte Carlo simulations for two value of exposure 
time of 3 200 s and 6 000 s, respectively. 



Element 


N lines 


<[X/H]> 


a 


N lines 


<[X/H]> 


a 






Exposure time of 3 200 s 




Exposure time of 6 000 s 


Nal 


0-1 


-3.170 ± 0.054 




0-1 


-3.178 ± 0.044 




Mgl 


0-1 


-2.527 ± 0.110 




0-1 


-2.550 ± 0.074 




Si I 


0-2 


-2.727 ± 0.182 


0.329 ± 0.125 


0-2 


-2.762 ± 0.102 


0.275 ± 0.072 


Cal 


6-16 


-2.602 ± 0.032 


0.061 ± 0.017 


7-15 


-2.621 ± 0.020 


0.043 ± 0.011 


Sell 


L4 


-3.085 ± 0.130 


0.105 ± 0.088 


1^1 


-3.091 ± 0.112 


0.106 ± 0.083 


Til 


5-17 


-2.559 ± 0.037 


0.076 ± 0.025 


5-16 


-2.587 ± 0.027 


0.048 ± 0.015 


Till 


10-22 


-2.564 ± 0.040 


0.094 ± 0.019 


10-23 


-2.590 ± 0.028 


0.067 ±0.014 


VI 


2-5 


-3.029 ± 0.062 


0.056 ± 0.047 


2-6 


-3.031 ± 0.039 


0.029 ± 0.026 


VII 


0-1 


-2.453 ± 0.093 




0-1 


-2.531 ± 0.067 




Cri 


2-4 


-3.011 ± 0.065 


0.046 ± 0.033 


24 


-3.015 ± 0.047 


0.034 ± 0.026 


Cm 


0-1 


-2.591 ± 0.198 




0-1 


-2.719 ± 0.127 




Mm 


2-10 


-3.025 ± 0.047 


0.068 ± 0.029 


2-9 


-3.036 ± 0.029 


0.044 ± 0.020 


Mnll 


0-1 


-3.067 ±0.118 




0-1 


-3.059 ±0.081 




Fel 


38-64 


-2.996 ±0.018 


0.078 ±0.011 


38-63 


-3.019 ±0.012 


0.053 ± 0.008 


Fell 


2-6 


-2.928 ±0.076 


0.101 ±0.045 


2-6 


-2.955 ±0.057 


0.071 ± 0.032 


Col 


1-8 


-3.044 ± 0.058 


0.071 ± 0.045 


1-7 


-3.047 ± 0.037 


0.047 ±0.275 


Nil 


1-5 


-2.930 ± 0.064 


0.048 ± 0.043 


1-5 


-2.972 ± 0.042 


0.029 ± 0.024 


Sril 


0-1 


-3.093 ±0.104 




0-1 


-3.091 ± 0.073 




Yll 


0-1 


-3.023 ±0.117 




0-1 


-3.019 ±0.082 




Zrll 


0-1 


-2.891 ±0.106 




0-1 


-2.939 ±0.071 




Ball 


1-2 


-3.099 ±0.107 


0.068 ± 0.059 


2 


-3.106 ±0.078 


0.051 ± 0.039 


La II 


04 


-2.792 ±0.170 


0.127 ±0.154 


04 


-2.874 ±0.122 


0.090 ±0.114 


Cell 


0-1 


-2.608 ±0.122 




0-1 


-2.651 ±0.091 




Ndll 


0-1 


-2.903 ±0.107 




0-1 


-2.964 ± 0.070 




Eull 


0-1 


-2.995 ± 0.084 




0-1 


-3.024 ± 0.062 





indicator of the error than the line-to-line scatter. The S/N 
ratio in the figures is derived in the centre of the blue CCD. 
The relations between the S/N ratios in the various posi- 
tions on the CCDs depends on the exact instrument design, 
and in this specific case we refer to the phase A 4MOST de- 
sign (Mignot et al. 2012). As expected, the line-to-line scat- 
ter decreases as the S/N ratio increases; the scatter on the 
average abundance and on the line-to-line scatter decreases 
also as the S/N increases. 

We then fit a second order polynomial on the line-to-line 
scatter as a function of S/N for those elements, for which 
the abundance ratios were derived from at least six lines 
(Si I, Cal, Til, Vl, Fel, Fell, Nil). The results are shown 
in Table[l~4l These values can be interpolated to reach the 
practical request that, in order for the line-to-line scatter to 
be smaller than 0.1 dex for all elements, a S/N > 50 per 
pixel at the centre of the blue CCD needs to be reached. 

For the cases where only a single line, or few lines, are 
detectable in the spectral range, the r.m.s. scatter on the av- 
erage value is a better measure for the abundance precision, 
as argued above. In analogy to above, we fitted a second 
order polynomial to the scatter on the average A(X), as de- 
rived from the Monte Carlo events, vs. S/N ratio. These 
results are shown in Table[l~4lfor all species. Hence, the con- 
straint S/N = 50 per pixel (again, in the centre of the blue 
CCD) for all species will yield ct^a(x)) < 0.1 dex, except 



for Cr II and La II for which we find a limiting precision of 
0.11 and 0.12 dex, respectively. 



0.15 - 



0.10 



0.05 - 



0.00 



t 1 1 1 1 1 1 1 1 1 1 r- 

o Fel 
x Cal 
■ Nil 



40 



60 

(S/N)/pixel 



80 



Fig. 4 The average line-to-line scatter for three exemplary 
elements with a large number of measurable lines, as a func- 
tion of the S/N ratio per pixel at the centre of the blue 
range. 
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Table 14 Coefficients of a 2nd order polynomial fit to the line-to-line scatter and <7 aV erage as a function of S/N ratio, for 
the case 5500/4.0/0.0. 



Element 


ao 


line-to-line scatter 

ai 


ai 






ao 








^"average 

ai 




02 








Nal 
















0.137 








-0.224 


X 


10" 


2 


0.135 


X 


10" 


4 


All 
















0.144 








-0.253 


X 


10" 


2 


0.146 


X 


10" 


4 


Mgl 
















0.908 


x 


10" 


1 


-0.153 


X 


10" 


2 


0.864 


X 


10" 


5 


Si I 


0.121 


-0.188 x 


10" 


2 


0.109 x 


10" 


4 


0.592 


x 


10" 


1 


-0.850 


X 


10" 


-3 


0.424 


X 


10" 


5 


Sill 
















0.190 








-0.341 


X 


10" 


-2 


0.207 


X 


10" 


4 


Cal 


0.993 x 10" 1 


-0.153 x 


10" 


2 


0.899 x 


10" 


5 


0.533 


X 


10" 


1 


-0.862 


X 


10" 


-3 


0.471 


X 


10" 


5 


Sell 
















0.133 








-0.212 


X 


10" 


2 


0.111 


X 


10" 


4 


Til 


0.144 


-0.243 x 


10" 


2 


0.151 x 


10" 


4 


0.619 


X 


10" 


1 


-0.103 


X 


10" 


2 


0.592 


X 


10" 


5 


Till 
















0.131 








-0.257 


X 


10" 


-2 


0.163 


X 


10" 


4 


VI 


0.902 x 10 _1 


-0.456 x 


10" 


3 


0.132 x 


10" 


-5 


0.576 


X 


10" 


1 


-0.106 


X 


10" 


2 


0.650 


X 


10" 


5 


Cri 
















0.142 








-0.264 


X 


10" 


-2 


0.159 


X 


10" 


4 


Cm 
















0.318 








-0.583 


X 


10" 


2 


0.338 


X 


10" 


4 


Mm 
















0.734 


X 


10" 


1 


-0.906 


X 


10" 


3 


0.410 


X 


10" 


5 


Fel 


0.125 


-0.156 x 


10" 


2 


0.934 x 


10" 


-5 


0.195 


X 


10" 


1 


-0.276 


X 


10" 


3 


0.151 


X 


10" 


5 


Fell 


0.187 


-0.248 x 


10" 


2 


0.155 x 


10" 




0.911 


X 


10" 


y 


-0.153 


X 


10" 


■> 


0.855 


X 


10" 




Col 
















0.118 








-0.225 


X 


10" 


-2 


0.138 


X 


10" 


4 


Nill 


0.126 


-0.207 x 


10" 


2 


0.117 x 


10" 


-4 


0.577 


X 


10" 


1 


-0.886 


X 


10" 


3 


0.438 


X 


10" 


5 


Yll 
















0.182 








-0.320 


X 


10" 


-2 


0.191 


X 


10" 


4 


Zrll 
















0.225 








-0.462 


X 


10" 


2 


0.297 


X 


10" 


4 


Ball 
















0.214 








-0.348 


X 


10" 


-2 


0.180 


X 


10" 


4 


La II 
















0.300 








-0.513 


X 


10" 


2 


0.288 


X 


10" 


4 


Nd II 
















0.221 








-0.365 


X 


10" 


-2 


0.210 


X 


10" 


4 


Eull 
















0.127 








-0.218 


X 


10" 


-2 


0.125 


X 


10" 


4 



0.15 



0.10 



A 
>< 



< 

V 



0.05 - 



0.00 



-i 1 1 1 1 1 1 1 1 1 1 — 

x Zrll 

■ ■ Ball 

o Eull 



40 



60 

(S/N)/pixel 



80 



Fig. 5 The scatter on the 1 000 Monte Carlo events on the 
average abundance for three representative elements with 
only few available lines, as a function of the S/N ratio per 
pixel at the centre of the blue range. 



5.5 Radial velocity precision 

We determined radial velocities by a cross-correlation of 
our noise-injected synthetic spectra, for different values of 
resolution and S/N ratio, against the synthetic spectra as 
templates. This study does not take into account any poor 
knowledge on atomic data (laboratory wavelengths and os- 
cillator strength), possible mismatches between the stel- 
lar spectrum and the cross-correlation mask, nor zero-point 



considerations. The mismatch between the observed spec- 
trum and the template will mostly depend on the type of 
star. Figure 1 of Nidever et al. (2002) illustrates this effect of 
the precision on radial velocity. For instance, measurements 
made with a Solar template lead to an mismatch, thus er- 
ror, that is increasing much faster towards F-type stars than 
towards K-type stars. The mismatch has also a greater ef- 
fect for hot stars (Verschueren, David, & Vrancken 1999). 
The zero point is directly dependent on the accuracy of the 
instrumental wavelength calibration and/or the manner in 
which it is corrected for using standard sources (stars, aster- 
oids, telluric rays, ...). 

In practice, we pursued three different approaches, the 
results of which are listed in Table[i~5land which are labelled 
as: (i) "Gauss": strong lines (Call-K and H; FFy and H<5) 
were masked out and the cross correlation function (CCF) 
peak was fitted with a Gaussian; (ii) "PolyM": strong lines 
are masked out as above, but the CCF peak was fitted with 
a third-order polynomial; (iii) "Poly": no line-masking has 
been applied and the CCF peak was determined by third- 
order polynomial fit as in (ii). 

We cannot see any tangible difference between the two 
spectral ranges we probed; thus we conclude that the Ca II- 
K and -H lines do not affect the radial velocity determina- 
tion in any way. The precision in radial velocity determined 
here is much higher than the one expected from Gaia (see 
Katz at al. 2004) for the entire sample of 4MOST stars. The 
reason is that the high resolution of 4MOST is larger than 
the one of Gaia. Furthermore, the blue range of the 4MOST 
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Table 15 Radial velocity errors. The two ranges are similar to our selected blue range but in one case Call-K is in, in 
the other not, to check the impact of this strong line in the radial velocity determination. 



T cff /logg//[Q/Fe]/[M/H] 
[K] / [cgs] / [kms" 1 ] / [dex] / [dex] 


R 


S/N 

r ■ 1— In 

[pixel J 


cry [km s 

(390-450 nm) 


(400-460 nm) 




Gauss 


PolyM 


Poly 


Gauss 


PolyM 


Poly 


4700 / 2.0 / 0.4 / 


-2.0 


15 000 


100 


0.022 


0.015 


0.014 


0.022 


0.015 


0.014 


4700 / 2.0 / 0.4 / 


-2.0 


15 000 


50 


0.044 


0.030 


0.027 


0.041 


0.029 


0.027 


4700 / 2.0 / 0.4 / 


-2.0 


15 000 


25 


0.084 


0.057 


0.052 


0.089 


0.062 


0.056 


4700 / 2.0 / 0.4 / 


-2.0 


20000 


100 


0.016 


0.011 


0.010 


0.016 


0.011 


0.011 


4700 / 2.0 / 0.4 / 


-2.0 


20000 


50 


0.031 


0.022 


0.021 


0.031 


0.022 


0.022 


4700 / 2.0 / 0.4 / 


-2.0 


20000 


25 


0.061 


0.044 


0.041 


0.062 


0.044 


0.041 


4700 / 2.0 / 0.4 / 


-2.0 


25 000 


100 


0.013 


0.010 


0.009 


0.013 


0.009 


0.009 


4700 / 2.0 / 0.4 / 


-2.0 


25 000 


50 


0.025 


0.019 


0.018 


0.025 


0.019 


0.018 


4700 / 2.0 / 0.4 / 


-2.0 


25 000 


25 


0.050 


0.037 


0.034 


0.048 


0.036 


0.035 


4750/2.0/0.0/ 


0.0 


15 000 


100 


0.013 


0.009 


0.009 


0.013 


0.009 


0.009 


4750 / 2.0 / 0.0 / 


0.0 


15 000 


50 


0.026 


0.019 


0.018 


0.026 


0.019 


0.017 


4750 / 2.0 / 0.0 / 


0.0 


15 000 


25 


0.052 


0.037 


0.034 


0.051 


0.036 


0.035 


4750 / 2.0 / 0.0 / 


0.0 


20000 


100 


0.009 


0.007 


0.007 


0.009 


0.007 


0.006 


4750 / 2.0 / 0.0 / 


0.0 


20000 


50 


0.019 


0.014 


0.013 


0.019 


0.014 


0.013 


4750 / 2.0 / 0.0 / 


0.0 


20000 


25 


0.037 


0.028 


0.025 


0.038 


0.027 


0.025 


4750 / 2.0 / 0.0 / 


0.0 


25 000 


100 


0.008 


0.006 


0.005 


0.008 


0.006 


0.005 


4750 / 2.0 / 0.0 / 


0.0 


25 000 


50 


0.015 


0.011 


0.011 


0.015 


0.011 


0.011 


4750/2.0/0.0/ 


0.0 


25 000 


25 


0.030 


0.023 


0.021 


0.031 


0.023 


0.021 



Table 10 Abundances results and uncertainties de- 
rived from 1 000 Monte-Carlo simulations for the 
(5500/4.0/-2.0) model and r = 14 mag. The expo- 
sure time, calculation method, and the meaning of the 
columns are the same as in Table|9] 



Element 


N lines 




<[X/H]) 




a 




Nal 


0-2 


-1 


.856 ±0.565 





.604 ± 


.425 


Mgl 


0-2 


-1 


.645 ± 0.077 





.077 ± 


.062 


Si I 


2-8 


-1 


.561 ±0.078 


0. 


.105 ±0 


.091 


Cal 


8-18 


-1 


.619 ±0.032 


0. 


.069 ± 


.016 


Sen 




-1 


.974 ± 0.302 





.247 ± 


.328 


Til 


2-11 


-1 


.561 ±0.062 





.108 ±0 


.049 


Till 


2-7 


-1 


.593 ± 0.080 


0. 


.081 ±0 


.041 


VI 


1-5 


-1 


.985 ±0.126 


0. 


.101 ±0 


.132 


Mm 


2-12 


-2 


.056 ± 0.063 


0. 


.107 ±0 


.038 


Fel 


40-71 


-2 


.006 ± 0.021 


0. 


.135 ±0 


.027 


Fen 


2-8 


-2 


.006 ± 0.078 


0. 


.102 ±0 


.052 


Co I 


1-7 


-2 


.068 ± 0.064 


0. 


.057 ±0 


.031 


Nil 


0-5 


-1 


.814 ±0.235 


0. 


.175 ±0 


.229 


Sril 


0-1 


-2 


.207 ±0.118 








Yll 


1-2 


-2 


.062 ± 0.082 


0. 


.047 ±0 


.045 


Zrll 


0-1 


-1 


.439 ± 0.081 








Ban 


1-2 


-2 


.038 ±0.111 


0. 


.057 ±0 


.050 


La II 


0-1 


-1 


.522 ±0.118 









Table 12 Abundances results and uncertainties derived 
from 1 000 Monte-Carlo simulations for the (4500/2.0/0.0) 
model. 

Element N lines ([X/H]) a 

Oi 
Nal 
Mgl 
All 
Si I 
Cal 
Scl 
Sc II 
Til 
Till 
VI 
Cri 
Mm 
Fel 
Fell 
Col 
Nil 
Yll 
Zrl 
Ball 
La II 
Prll 



1 


-0. 


.085 ± 0. 


.083 












0-2 


-0. 


.101 ±0. 


.043 





.050 


± 





.036 


1-2 


0. 


.032 ± 0. 


.028 





.017 


± 





.014 


2-2 


-0. 


.028 ± 0. 


.055 





.036 


± 





.027 


4-11 


-0. 


.071 ±0. 


.038 





.061 


± 





.019 


3-9 


-0 


.018 ±0. 


.027 





.039 


± 





.014 


0-1 


-0. 


.116 ±0. 


.054 












2-5 


-0. 


.017 ±0. 


.046 





.048 


± 





.024 


8-19 


-0. 


.089 ± 0. 


.025 





.056 


± 





.012 


1 


-0. 


.088 ± 0. 


.112 












5-15 


-0. 


.081 ±0. 


.018 





.033 


± 





.009 


1-2 


-0. 


.072 ± 0. 


.043 





.053 


± 





.039 


0-2 


-0 


.330 ± 0. 


.084 





.063 


± 





.043 


41-64 


-0. 


.065 ± 0. 


.015 





.064 


± 





.007 


2-4 


-0. 


.042 ± 0. 


.062 





.070 


± 





.037 


4-10 


-0. 


.042 ± 0. 


.022 





.040 


± 





.013 


7-15 


-0. 


.069 ± 0. 


.028 





.056 


± 





.014 


0-1 


-0 


.021 ± 0. 


.045 












2-4 


-0. 


.089 ± 0. 


.047 





.027 


± 





.020 


0-1 


-0 


.051 ±0. 


.083 












2-3 


-0. 


.026 ± 0. 


.066 





.028 


± 





.023 


1 


-0. 


.006 ± 


.107 













spectrograph contains a wealth of spectral lines, even in low 
metallicity stars, which facilitates the derivation of radial 

velocity at high precision. an order of magnitude. In practice, we expect constraints 

Bearing in mind that we are aiming at a precision of to rather come from other technical limitations of the spec- 

lkm s , these results suggest that the limitation to radial trograph (mechanical and thermal instabilities, photocentre 

velocity precision does not stem from S/N ratio, resolu- motion, scrambling achieved by the fibre, fibre cross-talk), 

tion or spectral coverage, as the errors in radial velocity in If such aspects can be kept under control, our results indi- 

our tests were always better than this specification by over cate that a precision in radial velocity of better than lkms -1 
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Table 13 Signal-to-noise ratio of the simulated spectra of 
r = 14 and exposure time of 3 600 s, at three wavelengths 
for both CCDs: 400, 420, and 450 nm for the blue and 590, 
620, and 670 nm for the red CCD. 



Spectrum 






S/N per pixel 








400 


420 


450 


590 


620 


670 


4500/2.0/0.0 


33 


37 


30 


70 


77 


50 


4500/2.0/-3.0 


29 


37 


28 


71 


79 


52 


5500/4.0/0.0 


44 


51 


36 


72 


78 


49 


5500/4.0/-2.0 


42 


50 


35 


72 


79 


49 



should be achievable. In all cases, the unmasked parabolic 
fits yield the smallest uncertainties. 

A consequence of this analysis is that, as long as we aim 
at a precision of the order of lkm s~\ the wavelength range 
and resolution of a high-resolution spectrograph should be 
rather driven by the optimisation of the chemical analysis, 
as elaborated in Section 5.3. 

6 Summary 

With accurate information for 1% of Galactic stars, Gaia 
will allow a breakthrough in our understanding of the MW. 
But this new knowledge would not be complete without the 
additional information that can be provided by high reso- 
lution spectroscopy. Anticipating the first data release of 
Gaia, and in the perspective of exploiting optimally our new 
insights on the Milky Way, there is a primary need for a 
follow-up ground-based spectroscopic survey. With this in 
mind, we explored the spectral quality such an instrument 
needs to deliver to achieve meaningful and accurate scien- 
tific results on stellar abundances. 

Based on extensive tests carried out on synthetic spectra 
over a broad range of stellar parameters we argue for the 
following, optimal and yet technically feasible design: 

- Wavelength ranges: 

blue arm: [395;459]nm, with eventually any gap 
between detectors avoiding to split the G-band , 
[429;432] nm, in the two detectors; 
red arm: [587;673] nm, with, if necessary, a gap around 
635 nm. 

- Resolving Power: 20000 on average across the above 
wavelength range, but never smaller than 18 000. 

- Pixel-Size is not a strong limitation, and we sug- 
gest to have an optimal sampling with two elements. 
As good compromises we envisage as values for Re- 
solving Power/Sampling: 25 000/2.0-2.5, 20000/2.5, 
15 000/2.5-3.0. 

An instrument with such a design will allow to derive 
abundances for the elements listed in TableQ] It can provide 
a precision in the abundance determination of about 0. 1 dex, 
for the majority of the elements. Tables IUEI and [6] can be 
used to estimate the typical accuracies from the single lines 
as a function of line strength and S/N ratio. 



None of the above has a major influence on the reach- 
able precision on radial velocity. Moreover, this is limited 
by systematics in the instrument design such as mechanical 
stability, which is beyond the scope of the present spectral 
feasibility study. 
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2i/2Tr72 2.355 
If we introduce the Sampling, S, as S = p^™™ we obtain 



N ; 



In y /In 2 w 1.2 x sJlrij. 



The number of pixels that define the line is proportional to the 
square root of the logarithm of the strength of the line itself. 

For a sampling of 2.0, a strong line with a residual intensity of 
0.1 is thus defined by only ~5 pixels, while a weak line of resid- 
ual intensity 0.9 covers less than 4 pixels. With the 2.5 sampling 
of the centre of the ranges, strong and weak lines are defined by 
more than 6 and less than 5 pixels, respectively. Likewise, at a 
sampling of 3.0 these numbers become less than 8 and more than 
5, respectively. This should be considered a lower limit, however, 
since strong lines are better described with a wider Voigt profile. 



A Relation of pixel vs. line strength 

In what follows, we represent an absorption line as a Gaussian 
profile: 

(A-A ) 2 

F(X) = 1 - Ae 2.2 . 

We take the relative depression of the line core A as a measure 
of the line strength. We consider a line as well-defined when it 
covers a certain number of pixels from the core to the continuum. 
In particular, we parameterise the extent of the line by chosing a 
suitable (small) number 8 such that the feature transits into the 
continuum if F > 1 — 5. Obviously, 8 is dependent on the S/N 
level, but it is typically of the order of 5 —0.01 for continuum 
normalized spectra. 

Here we are interested in determining a relation between the 
line strength, A, and the number of pixels it covers. From the above 
relation 



it follows that 



(A-A )^ 

Ae 2.2 < i 



(A — A ) 2 < 2a 2 In 4, 
o 



which results in a full width of the line of 



2|A, - A ; = 2 , 2\nj. 

where Ai is the wavelength where the line profiles transits into the 
continuum. 

The number of pixel defining the line profile is thus: 



N = 2 



PixelSize 



Now we intend to make a connection between a and spectral 
resolving power given by the FWHM of the spectrograph's instru- 
mental profile AAfwhm • At the resolution of 4MOST, stellar lines 
are typically not resolved, and the observed line profile largely re- 
flects the instrumental profile so that 
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